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ABSTRACT 

We use a high-resolution cosmological simulation to predict the distribution of H I Lya emission from the low- 
redshift (z < 0.5) intergalactic medium (IGM). Our simulation can be used to reliably compute the emission 
from optically thin regions of the IGM but not that of self-shielded gas. We therefore consider several models 
that bracket the expected emission from self-shielded regions. Most galaxies are surrounded by extended 
(> 10 2 kpc) "coronae" of optically thin gas with Lya surface brightness close to the expected background. Most 
of these regions contain smaller cores of dense, cool gas. Unless self-shielded gas is able to cool to T < 10 4A K, 
these cores are much brighter than the background. The Lya coronae represent "cooling flows" of IGM gas 
accreting onto galaxies. We also estimate the number of Lya photons produced through the reprocessing 
of stellar ionizing radiation in the interstellar medium of galaxies; while this mechanism is responsible for 
the brightest Lya emission, it occurs on small physical scales and can be separated using high-resolution 
observations. In all cases, we find that Lya emitters are numerous (with a space density ~ 0.1 Mpc~ 3 ) and 
closely trace the filamentary structure of the IGM, providing a new way to map gas inside the cosmic web. 
Subject headings: cosmology: theory - galaxies: formation - intergalactic medium - diffuse radiation 



1. INTRODUCTION 

One of the most striking images to emerge from cosmo- 
logical simulations is the "cosmic web". Matter collapses 
into a web of moderately overdense sheets and filaments, with 
galaxies and galaxy clusters forming through continued col- 
lapse at their intersections. This picture explains both the 
qualitative distribution o f galaxies in redshift surveys (e.g., 
de Lappare nt et al.lll986|) and many of the characteristics of 
the Lya forest (see e.g. Rauch 1998). However, mapping the 
cosmic web remains a difficult venture. Lya forest spectra 
probe only one-dimensional lines of sight, while galaxy red- 
shift surveys find only the small fraction of baryons embedded 
inside galaxies, offering an indirect picture of the gaseous fil- 
aments. 

In this Letter, we use a cosmological simulation to show 
that surveys of hydrogen Lya (A 1 2 1 6 A) emission can provide 
a powerful method to map the cosmic web. Existing analytic 
estimates suggest that the surface brightness of de nse portions 
of the IGM can be substantial (Hogan & Wevmann 1987; 
iGould & Weinbergl ll996). Such regions generally lie within 
filaments but outside of galaxies, so Lya emission offers a 
more direct map of the gas distribution in filaments than do 
galaxy surveys. In fact, Lya emission from galaxi es has been 
used to detect a filament at high redshift ( Mo ller & Fvnbol 
2001); here we focus on intergalactic emission from low red- 
shifts, where the distances and sky background are relatively 
small and where Lya maps of the cosmic web can be com- 
pared to existing galaxy surveys. We also show that Lya- 
emitting gas constitutes an intrinsically interestin g phase: gas 
that is cooling aro und collapsed halos (Haiman et al.I l2000t 
iFardal et alJl200l() . The distribution of Lya emission there- 
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fore provides a window onto the growth of bound objects. 

2. SIMULATION AND ANALYSIS 

We perfo rm our analysis us ing the G5 cosmological simula- 
tion of Springel & Hernquist (2003b). This smoothed-particle 
hydrodynamics simulation included a multiphase description 
of star formatio n which incorporates a p rescription for galac- 
tic winds (Springel & Hernquist 2003a). The simulation as- 
sumed a ACDM cosmology con sistent with the m ost recent 
cosmological observations (e.g., Sperg eTet al.ll2003T) . It has 
a box size of 100/T 1 comoving Mpc (where the Hubble con- 
stant Ho = 100/i km s _l Mpc" 1 and h = 0.7 in the simulation), 
a spatial resolution of 8/i _1 comoving kpc, and a baryonic par- 
ti cle mass of 3.26 x IP 8 /? " 1 M . 

Furl anetto et al.l d2003l hereafter F03), who focus on UV 
emission from metals in the IGM, describe the main compo- 
nents of our analysis procedure; we briefly summarize them 
here. We first compute two grids of the Lya emissivity e a 
as a function of hydrogen de nsity riw and ga s temperature T 
using CLOUDY 96 (beta 5, Ferland 2003), assuming ion- 
ization equilibrium. F or the first grid (£«'), we apply the 
lHaardt & MadaulJ20 01 ) ionizing background, which includes 
radiation from galaxies and quasars (our main results are in- 
sensitive to the background choice; F03). For the second 
grid (e c ° 11 ), we neglect all photoionization processes; we use 
this grid to calculate the emissivityof self-shielded gas. 5 We 
compare the two grids in Figure ^ The solid curve shows 
e c ° n (note that e a oc njj when collisional processes dominate), 

and the other curves show e P d for several different densities. 
We see that e c ° n drops rapidly when the gas recombines at 

T w 10 4 ' 2 K. In contrast, e p d is still large at low temperatures 
because the gas remains photoionized. As the density and 
temperature increase, collisional processes become dominant 

and et -> e c Q o11 . 
We wish to compute the emission from a slice of the sim- 

5 Note that the electron density (and hence e™ ) becomes sensitive to the 
metallicity for T < 10 4 K. However, this has negligible effects on our results 
because at these temperatures the emissivity is vanishingly small. 
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FIG. 1. — The normalized Lya emissivity e a /n^, in units of erg cm 3 s . 
The solid curve includes only collisional processes. The dashed, dot-dashed, 
and dotted curves allow photoionization and show the emissivity for gas with 
«H = 10"*, 10~ 4 , and 10~ 2 cm -3 , respectively. 



ulation with fixed z, Az, and angular size. We first ran- 
domly choose the volume corresponding to this slice (see 
F03). We classify particles inside the slice as "optically 
thin" or "self-shielded." The latter set includes those dense 
and cool particles with a large enough neutral fraction to 
shield themselves from the ionizing background. We approx- 
imate self-shielding by assuming that all gas with «h > «ss 
and T < r ss is optically thick. We set fiducial thresholds at 
« ss = 10~ 3 cm" 3 and r ss = 10 4 5 K; radiative transfer calcula- 
tions f or self-gravitating clouds like those described in Schave 
fcOOll) find that self-shielding becomes important at about this 
densi ty, assuming that it remains cool (see also iKatz et all 
As shown in F03, most particles in this regime lie on a 
well-defined curve in the n^-T plane, beginning at moderate 
overdensity and T ~ 10 4 25 K and asymptotically approach- 
ing 10 4 K as the density increases to «h ~ 0.1 cm" 3 . This 
locus represents gas that is cooling and falling into the centers 
of halos. A small but important fraction of particles (those 
that have recently been heated by shocks) have slightly higher 
temperatures. Once we have made this classification, we as- 
sign emissivities to optically thin particles using the e P d grid 
and (in our fiducial model) to self-shielded particles using the 
e™ 11 grid. 6 Unfortunately, the emissivity of self-shielded par- 
ticles is rather uncertain because the simulation does not in- 
clude radiative transfer, metal line cooling, or local ionizing 
sources and because we assume ionization equilibrium . Fo r 
this reason, we examine their emission more closely in 33.21 
Active star formation also produces substantial Lya emis- 
sion when gas in the host galaxy absorbs ionizing radiation 
from stars and subsequently recombines. In the simulation, 
star formation occurs in any gas particle whos e density ex- 
ceeds «h = 0.129 cm -3 (jSpringel & Hernauist 2003a). For 
each such particle, we convert the star formation rate (SFR) 
reported by the simulation to a Lya luminosity via L a = 
10 42 (SFR/M©yr _1 ) erg s" 1 , which is accurate to within a fac- 
tor of a few for meta llicities between 0.05 Z Q < Z < 2 Z Q 
jLeitherer et all 1999). assuming tha t ~ 2/3 of ionizin g pho- 
tons are converted to Lya photons (Osterbrock 1989). The 
assumption of case-B recombination should be valid provided 
that most ionizing photons are absorbed in the dense inter- 
stellar medium of the galaxy. The actual Lya luminosity 

6 Note that this is conservative in the sense that we neglect photoionization 
by hard UV/X-ray background photons to which the gas is still optically thin. 



of a galaxy depends strongly on the distribution of ionizing 
sources, the escape fraction of ionizing photons, the pres- 
ence of dust, and the kinematic structure of the gas (e.g., 
iKunth et"afll2003l) . so this should be taken as no more than 
a representative estimate. Note that any star formation un- 
resolved by the simulation could significantly change our re- 
sults. However, Sprinael & Hernauist (2003b) show that the 
global SFR in the simulation matches observational estimates 
well, so such difficulties are probably not severe. 

Finally, we locate each particle on a pixelized map, convert 
its emissivity to surface brightness, and smooth the result with 
a Gaussian filter of full width at half maximum Ad (F03). 

3. RESULTS 

Figure|2jj shows a thin slice (Az = 10~ 3 ) of our simulated 
universe at z = 0.15. The panel shows a large volume with 
coarse resolution; the colorscale includes the full dynamic 
range of surface brightness $. Note how filaments stand 
out clearly in the map. Figure 2b shows part of a filament 
with a finer resolution (~ 13 h~ l physical kpc). We include 
only those pixels with $ > 10 photons cm s" 1 sr" 1 ; this is 
a fact or ~ 10 smaller t han some observational l lBrown et alJ 
2000) and theoretical (Haardt & Madaul l2"00"ll) estimates of 
the diffuse background at 1200-2000 A. While the mean 
emission from filaments lies well below the background, we 
find that numerous discrete regions trace their course: the 
space density of these Lya emitters is ~ 0.1 h 3 Mpc~ 3 . Wide- 
field Lya observations are thus a promising technique to map 
out the cosmic web. Note that decreasing the threshold to 
$ > 1 photons cm" 2 s" 1 sr" 1 enlarges the coronae but does not 
significantly increase their number density. 

Lya emission occurs over a wide range of physical scales: 
generally, emitting regions have a central, high surface bright- 
ness (i> > 10 3 photons cm" 2 s" 1 sr" 1 ) core of size ~ 30- 
50/T 1 kpc surrounded by an irregular, low surface brightness 
"corona" that can extend to scales > 10 2 kpc. Thus, the Lya 
emission is much more extended than the galaxies, indicating 
that this view of the cosmic web traces truly intergalactic gas. 

Figure|3]shows histograms of the pixel flux probability dis- 
tribution function (PDF) dP/dlog$. The solid curve in each 
panel shows our fiducial model (using e c ° u for self-shielded 
gas). Figure 3a shows the PDF to very small $; we see a broad 
peak centered at $ ~ 10~ 3 photons cm" 2 s" 1 sr" 1 and a tail ex- 
tending to much higher surface brightness. As described in 
Paper I, the location of the peak is determined simply by ion- 
ization equilibrium i n the mean-density IGM and can be es- 
timated analytically (Gould & Weinberg 1996). Note that the 
PDF depends strongly on A9, because the emitting regions 
have structure on small scales. 

Figure 3b shows that the PDF evolves very little from z = 
0.07 to z = 0.45. The fraction of bright pixels increases slowly 
with redshift because both the density and temperature of the 
IGM increase with redshift. 

3.1. Lya Photons From Star Formation 

Figure 2d shows the Lya emission from star formation. 
While all star-forming regions are surrounded by intergalac- 
tic Lya emission, the converse is not true. We do find that 
Lya emission nearly always surrounds star particles, but of- 
ten these are dwarf galaxies or lack ongoing star formation. 
The exceptions are galaxies inside hot clusters, which do not 
always host Lya emission. Note that Lya emission produced 
through star formation is much more compact than the coro- 
nae (partly because we assume that all the ionizing photons 
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FIG. 2. — Maps of Lya surface brightness $ for z = 0.15 and Az = 10 3 (AA = 1.2 A). Panel (a) assumes an angular resolution of ~ 29"; the rest show the 
region outlined in green with 7.2" resolution (or ~ 13 A -1 kpc). Except for panel (a), we exclude pixels with $ < 10 photons cm~ 2 s" 1 sr~' . (a), (b): Fiducial 
model, (c): e a = for self-shielded gas. (d): Estimated Lya emission from star formation in the slice (such emission was not included in the other panels). 




FIG. 3. — The pixel surface brightness PDF in maps with z, = 0.15 (unless otherwise specified), Az = 10~ 3 (AA =1.2 A), and 13/j~' kpc resolution. In all panels, 
the solid curve shows our fiducial model, (a): e a = for self-shielded gas (dashed curve) and no self-shielding cut (dotted curve), (b): z = 0.07, 0.15, 0.3, and 
0.45 (dotted, solid, long-dashed, and short-dashed curves, respectively). Also shown is the Lya emission from star-forming gas (dot-dashed curve) at z = 0.15. 
(c): e a =0 for self-shielded gas, with n ss = 10" 4 , 10~ 3 , and 10~ 2 cm -3 (dotted, short-dashed, and long-dashed curves, respectively). 



are absorbed locally). Figure 3b compares the mechanisms in 
a more quantitative form. The emission from star-forming gas 
in the simulation (dot-dashed curve) dominates the brightest 
pixels but is unimportant for $ < 10 5 photons cm" 2 s" 1 sr" 1 . 
If a large fraction of the ionizing photons escape their host 
galaxy, the number of pixels affected by local star-formation 
will increase, but their mean brightness will also decrease. 
This could be particularly important near active galactic nu- 
clei (AGN), which have large ionizing fluxes and large escape 
fractions. AGN could "light up " large volumes aroun d their 
host galaxies in Lya emission (Haiman&Rees 2001). We 
note that the SFR, and hence Lya emission from s tar-forming 
gas, increases with redshift (Hernquist & Snring ell2003l) . 

3.2. Self-Shielded Gas 

The largest uncertainty in our calculation is probably the 
treatment of self-shielded gas. One problem is identifying 
which particles are shielded from the ionizing background; 
we find, however, that our results do not change significantly 
in the range 10" 4 cm" 3 < n ss < 10~ 2 cm" 3 , because most gas of 
these densities has T « 10 4 1 K, where eS « e c ° u (see Fig.Q. 
A more important concern is that we may have overestimated 



the temperature of self-shielded particles because the simula- 
tion erroneously allows the ionizing background to heat such 
particles and underestimates the cooling rates in optically 
thick regions. Radiative transfer calculations with CLOUDY 
suggest that the simulation may have overestimated the tem- 
perature of particles with «h ~ 10~ 3 cm" 3 by > 25%, which 
can be significant given the steep dependence seen in Fig- 
ure ^ The error increases as the gas density and metallic - 
ity increase. In a worst case scenario, all self-shielded gas 
would have T < 10 4 1 K, producing essentially no emission. 
In Figure 2c, we approximate this case by setting e a = for 
all self-shielded gas. Comparing to the fiducial model, we 
find that the core brightness declines dramatically, but the ex- 
tended coronae are unaffected. We conclude that the large- 
scale emission is due to optically thin gas while the brightest 
regions consist of gas at or near to the self-shielding thresh- 
old. The resulting PDF is shown as the short-dashed curve 
in Figure 3a; it essentially eliminates the bright tail of the 
distribution. Figure 3c shows that our results do depend on 
n ss if self-shielded gas has e a = 0. The PDF changes little 
for 10" 4 cm" 3 < n ss < 10" 3 cm , because such gas has rela- 
tively low density and hence small emissivity. However, it ap- 
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proaches the fiducial model if we set n ss = 10~ 2 cm -3 , because 
gas with «h > 10~ 2 cm" 3 has T « 10 4 K in the simulation 
and hence a small e°° . 

On the other hand, the simulation also neglects local ioniz- 
ing sources. We have already seen that most self-shielded par- 
ticles reside near galaxies; we would therefore expect the ra- 
diation field around these particles to exceed the background. 
In this case the emissivity is limited by the ionizing flux reach- 
ing the gas. We have esti mated the number of ionizing pho- 
tons from young stars in 33.11 although we did not attempt 
to model the spatial distribution of the resulting emission. In 
addition, dense gas may no longer be self-shielded, so the dif- 
fuse background will also contribute to the Lya emission. We 
estimate this extra component by assuming that no particles 
are self-shielded (dotted curve in Figure 3a). The number of 
bright pixels increases because eg > e™ 11 for T < 10 4 1 K. In 
either case, the net effect of local ionizing sources is to in- 
crease the Lya emission. 

4. DISCUSSION 

We study Lya emission from gas in the IGM using a high- 
resolution cosmological simulation. We find that galaxies 
are typically surrounded by "coronae" of Lya emission with 
bright, central cores (of diameter < 40 h kpc) surrounded by 
much larger (> lQ 2 h~ l kpc) regions with surface brightness 
near the background. While the number of Lya photons due 
to star formation in the host galaxies can exceed that emitted 
by the surrounding gas, the spatial extent of the star formation 
is typically much smaller. High angular-resolution observa- 
tions (probing physical scales < 15 h~ l kpc) should be able to 
separate the two components. Detection of the bright cores 
is feasible with existing technology, although the wide-field 
ultraviolet spectrographs most useful for these studies have 
not yet been built. The Galaxy Evolution Explorer (GALEX), 
with a large field of view but relatively low spectral resolu- 
tion (~ 10 A), may detect the brightest cores. Furthermore, 
the strong correlation between Lya emission and galaxies 
may enable a statistical detection of the signal in deep ex- 
posures. A project underway to construct a balloon-borne 
wide-field UV spectrograph will have a limiting sensitivity 
<£> ~ 500 photons cm" 2 s" 1 sr" 1 in a single night's observation 
(D. Schiminovich, private communication). 

We emphasize that the emission from self-shielded gas 
is uncertain, because the simulation does not include all of 
the relevant physics. We have conside red a range of mod- 
els bracketing these uncertainties (see 33.2t : in our fiducial 
model, dense, cool gas produces Lya photons only through 
collisional processes, but we also consider cases in which it 
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has zero emissivity or in which we allow photoionization. The 
choice has little effect on the large-scale, low surface bright- 
ness emission, which comes primarily from optically thin gas, 
but it does strongly affect the luminous cores., where most of 
the gas is shielded from the ionizing background. 

We have also neglected dust, which can efficiently destroy 
Lya photons. Because the emitting region is often consider- 
ably larger than the associated galaxy, much of the emitting 
gas may be relatively pristine. In the simulation, we find that 
the majority of particles with large Lya emissivity contain 
no metals (or dust), though note that the simulation does not 
include mixing between simulation particles. At higher red- 
shifts, several large- scale Lya-emitting "blobs" have already 
been observed (e.g., Steidel et al. 2000), limiting the amount 
of dust in these environments. 

The structure and relatively large spatial size of Lya emit- 
ters, together with their locations, suggest that they are "cool- 
ing flows" of gas onto halos. Detailed, high resolution Lya 
observations therefore offer a means to study the accretion 
and cooling of gas onto galaxies. Galactic winds will also af- 
fect the distribution of Lya emission. However, we find that 
the winds in our simulations have only minor effects on the 
statistics of Lya emitters. 

On larger scales, the Lya emitters are distributed along 
sheets and filaments. Because these objects are numerous 
(with a space density ~ 0.1/i 3 Mpc" 3 ), Lya line observations 
offer a new and powerful way to map the cosmic web. Figure 
|2] shows that such observations do not require high resolu- 
tion; wide-field observations could efficiently locate (in three 
dimensions) a large number of Lya emitters. In contrast to 
galaxy surveys, this approach reveals the location of baryons 
outside of galaxies. Because Lya emission traces the cool gas 
accreting onto galaxies, it also selects a different population of 
dark matter halos, including many with few (or even no) asso- 
ciated stars. Follow-up observations at high resolution would 
allow us to compare Lya emitters to galaxy surveys and learn 
how galaxies accrete gas at the present day. 
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